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Abstract. Based on results from femtosecond time-resolved photoemission, we 
compare three different methods for determination of the electron-phonon coupling 
constant A in Eu and Ba-based 122 FeAs compounds. We find good agreement between 
l/"~j all three methods, which reveal a small A < 0.2. This makes simple electron-phonon 

m mediated superconductivity unlikely in these compounds. 
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1. Introduction 

Albeit an enormous amount of research has been conducted in the past few years 
on the Fe pnictide high-temperature superconductors (HTSCs) [H [2], the search for 
the superconducting pairing mechanism in these materials is still ongoing. Besides 
other excitations such as spin fluctuations [3], the electron-phonon (e-ph) coupling 
responsible for Cooper pairing in conventional superconductors in the Bardeen-Cooper- 
Schrieffer (BCS) theory is considered a potential candidate. Therefore, the quantitative 
determination of the e-ph coupling constant A in the Fe pnictides is of particular 
interest. Several methods to determine A have been established in the literature and 
successfully used in various studies. In the energy domain, angle-resolved photoemission 
spectroscopy (ARPES) analyzes the renormalization of the single-particle spectral 
function by the electronic self energy S, which includes the e-ph coupling. However, 
such type of analysis requires very high sample and data quality, which is difficult to 
obtain in the Fe pnictides, and only few ARPES studies of FeAs compounds were able 
to determine effects of S [H [5j [6]. In addition, in thermal equilibrium, contributions 
from other degrees of freedom to E are often difficult to disentangle from e-ph coupling. 

In contrast, femtosecond (fs) time- resolved spectroscopies allow to separate e-ph 
scattering from other relaxation channels like e.g. electron-electron (e-e) scattering or 
heat diffusion due to their different intrinsic timescales out of thermal equilibrium [TJ 
[HJ [9]. Here, time-resolved ARPES (trARPES) is used, which combines the energy and 
momentum resolution of ARPES with fs time resolution into a powerful tool to directly 
investigate the dynamics of the electronic structure in a non-equilibrium state [TQl ITT] . 

The e-ph relaxation in metals after photoexcitation has been successfully described 
using the two-temperature model (2TM) [TJ El El H21 [131 El- This model describes 
a system of two coupled heat baths for the conduction electrons and the ion lattice 
with temperatures T e and TJ by coupled rate equations for T e and TJ, respectively. One 
key assumption of this model is that e-e and phonon-phonon (ph-ph) scattering occurs 
on a much faster timescale than the e-ph scattering, leading to a thermal distribution 
within each subsystem. However, in correlated electron systems like the cuprate or Fe 
pnictide HTSCs, e-ph scattering might occur on similar timescales [IS], which poses 
some questions on the applicability of the 2TM for these materials and care has to 
be taken. Therefore, we analyze our data in terms of a suitable version of the 2TM 
and compare the results with two complementary methods for the determination of the 
e-ph coupling strength. We find that all three methods to analyze the e-ph coupling 
strength from time-resolved photoemission data yield results for the second moment of 
the Eliashberg coupling function A (u 2 ) which agree well within error bars and reproduce 
the trends for three different 122 FeAs compounds consistently. 
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Figure 1. (a) Scheme of the pump-probe experiment, (b) Time-resolved ARPES 
spectra of EuFe2As2 at T = 100 K and at normal emission for various pump- 
probe delays on a logarithmic intensity scale using an incident pumping fluence of 
F = 0.8mJ/cm 2 . Lines are fits to Fermi-Dirac distribution functions, (c) Electronic 
temperature T e and a fit to the 3TM. Temperatures of the hot phonons (T p ) and the 
rest of the lattice modes (T/) are shown as dashed and dash-dotted lines, respectively. 



2. Experiments and results 

For trARPES experiments, single crystals of EuFe2As2 and BaFe2As2 parent compounds 
and optimally doped BaFei.8sCoo.i5As2 (T c = 23K) were cleaved in ultrahigh vacuum 
(p < lCr 10 mbar) at T = 100 K where most measurements were carried out. The 
experimental setup is sketched in fig. Ilia): The output of a commercial regenerative 
Ti: Sapphire amplifier (Coherent RegA 9050) delivering ultrashort laser pulses at hv\ = 
1.5 eV photon energy with a pulse duration of 55 fs, and operating at 300 kHz is used 
to optically excite the samples (pump pulse). A time-delayed frequency-quadrupled 
probe pulse at hv-i = 6.0 eV photon energy with a pulse duration of 80 fs leads to the 
emission of photoelectrons, which are detected using an electron time-of-flight (TOF) 
spectrometer with an acceptance angle of ±3°. The energy resolution of 50meV is 
mainly determined by the spectral width of the probe pulses, and the overall temporal 
resolution is < 100 fs. For details see [9]. 
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2.1. Three- Temperature model 

Exemplary trARPES data of EuFe2As2 are shown in fig. [lib) for various pump-probe 
delays on a logarithmic intensity scale. After excitation, a pronounced distribution 
of excited charge carriers is formed at E — Ep > 0.1 eV, which deviates from the 
Fermi-Dirac distribution of the electronic system before excitation and which originates 
from hot, non-thermalized electrons. Subsequently, e-e scattering quickly leads to 
thermalization of these non-thermal electrons and to the formation of a hot thermalized 
electron distribution [9j. At t > 250 fs, the non-thermal contribution has decayed to less 
than 1% of the electron population. 

The temperature of the thermalized part of the electronic system T e can be 
extracted from the trARPES data by fitting a Fermi-Dirac distribution function to the 
high-energy cutoff of the transient spectra, multiplied by a phenomenological density- 
of-states function and convoluted with an instrumental resolution function [9j [H], as 
shown in fig. [lib) (solid lines). T e determined by the fitting is shown in fig. HTc) as a 
function of pump-probe delay. After the steep rise of the electronic temperature at zero 
pump-probe delay, we find a relaxation on two distinct timescales r a and Tp. Such a 
behavior indicates the selective coupling of hot electrons to a subset of strongly coupled 
phonon modes on the faster timescale r a and the subsequent energy transfer to the rest 
of the phonons by anharmonic ph-ph scattering, determined by Tp. Thus, the system is 
described by an extended version of the 2TM, which had been developed for the cuprate 
HTSCs [14|, where a similar relaxation dynamics was found, and which has been also 
used recently for the Fe pnictides pEElEE]. This three-temperature model (3TM) consists 
of a system of three coupled differential equations describing the temperature of the 
electrons, T e , of a hot phonon distribution T p and of the rest of the lattice modes, T\. 

f = -H<?.,T,) + §- (1) 



-f= +^H(T e ,T p )-^- 

dt C p Tp 



+ -^H(T e ,T p ) - ^ l - (2) 



d-/ / Gr> In 1] /„N 

d^ = + a^ ■ (3) 

The source term S describes the optical excitation, C e , C p and Ci are the specific heat 
capacities of electrons, strongly and weakly coupled phonon modes, respectively. The 
anharmonic decay of phonons is described by Tp. For the energy transfer from the 
electrons to the more strongly coupled phonons, the formula derived by Allen [7] is 
used: 

H{T n T,) = 7T (T, -T,)= m lf ] T < ~ T ' , (4) 

irk B J-e 

where k B is the Boltzmann constant. This relation allows determination of the second 
moment of the Eliashberg e-ph coupling function A (oj 2 ). 

A fit of a numerical solution of the 3TM to T e is shown in fig. [Tic) and yields a 
good agreement to the experimental data. Due to an ambiguity of the model parameters 
such as the electron and lattice specific heat capacities reported in the literature, a 
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Figure 2. (a) e-ph and e-e contributions to the electron decay rate I\ The e-ph 
contribution T e _ p h calculated in the Debye model with huju = 20meV and A = 0.3 at 
T = OK and T = 100 K increases up to ~ u>d and is constant above. The electronic 
contribution T e _ e calculated for j3 = 0.1 eV" 1 exceeds T e _ p f l only at higher energies. 
The shaded area marks the energy where e-ph scattering dominates, (b) Electron mean 
excess energy extracted from trARPES data at T = 100 K near kp within the energy 
window 0.07 eV < e < 0.3 eV marked in (a). Solid lines are fits to equation §6b. 



range of values for A (u 2 ) is retrieved, where the fits show similar good agreement to 
the data. Details of this analysis can be found in [17]; here, only the main results are 
discussed. From the 3TM, we find values of A (u> 2 ) = 56 — 65meV 2 for EuFe2As2, while 
in BaFei.85Coo.i5As2 it is slightly smaller (A (u 2 ) = 46 — 55meV 2 ) and even smaller in 
undoped BaFe 2 As 2 (A (u 2 ) = 30 - 46meV 2 ). 



2.2. Electronic excess energy 

The analysis within the 3TM was based on the energy relaxation within the thermalized 
part of the transient electronic distribution function and thus neglected the non-thermal 
electrons present at early delay times (see fig. pTb)). However, the analysis of the energy 
relaxation of these non-thermal, excited electrons can also provide information about 
the strength of e-ph coupling. 

The scattering rate of electrons excited at energy e = E—E-p above the Fermi level is 
in the self-energy formulation of many-body theory determined by F = hr^ 1 = 2ImE(e), 
where ImE is the imaginary part of the electronic self energy [TS]. Important energy- 
dependent contributions to F arise from e-e and e-ph scattering, where e-e scattering 
is considered to follow the quadratic energy scaling of the Fermi liquid theory, F e _ e = 
2/3 [(nksT) 2 + e 2 ] (fig. p2[a)). As the proportionality coefficient (3 is rather small in 
typical metals (in the order of 10 _2 eV _1 to 10 _1 eV _1 ) [THl [20], this contribution is 
negligible compared to the e-ph scattering at low enough excitation energies. The latter 
increases up to the maximal phonon energy hu} max and is constant above (fig. |2](a)). For 
e > hu max and T = K, F e _ ph results for an Einstein mode u to [2"T] 



e—ph 



TthXuJn 



(5) 
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Within the energy window between hu max and the crossover regime, where r e _ e becomes 
dominant (shaded area in fig. [21(a)), the rate of energy dissipation of an electron due to 
the emission of a phonon with energy Jvjjq is given by 

— = — = nh\uj , (6) 



which leads to a linear relaxation of the electron energy 

The rate of energy relaxation of excited electrons can be extracted from the 
experimental trARPES intensity I{e,t) by analyzing the mean excess energy, 

<A£(t)> = % \ . (7) 

However, this integral represents the mean excess energy within a selected partition 
of a transient distribution function and not the energy relaxation of single individual 
electrons. Thus, special care has to be taken with the determination of the integration 
boundaries, €q and e\. A careful investigation of the influence of the integration 
boundaries and the excitation fluence reveals that a reasonable choice of eo = 70 meV 
and ei = 300 meV allows to give a good estimate of the e-ph coupling strength, as 
detailed in [23J. 

Experimental data of EUES2AS2, BaFei.ssCoo.isAsa and BaFe2As2 near kp and at 
T = 100 K are shown in fig. [21(b) . To minimize lattice heating and the influence of 
the hot thermalized electron distribution, low excitation fluences of F ~ 50/iJcm -2 
have been used. The linear fits to equation (J6J) within the first 100 fs yield values of 
\(co 2 ) = 65(5) meV 2 for EuFe 2 As 2 , 50(3) me V 2 for BaFei. 85 Co .i 5 As2 and 34(6) me V 2 
for BaFe2As2- Here, error bars represent the numerical uncertainties of the fits, while 
the overall errors have to be considered larger due to the uncertainty of the integration 
boundaries discussed above and of the elevated lattice temperature. Nevertheless, 
we find a good agreement with the values obtained in the 3TM simulations and a 
considerably higher e-ph coupling in EuFe2As2 than in BaFei.85Coo.i5As2 and even lower 
A (u 2 ) in BaFe 2 As 2 . 

2.3. Temperature- dependent hole relaxation rates 

Finally, another estimate of the e-ph coupling strength can be gained from the 
temperature dependence of quasiparticle (QP) relaxation rates [3, [15], EH ES] . A recent 
theoretical investigation of the Boltzmann equation for e-e and e-ph interaction found 
an analytic solution for the temperature dependence of the QP relaxation rate r [24J . 
In the limit of a bad metal at elevated temperatures, r depends linearly on the lattice 

temperature T [TBI 1241 [25] : 

2irk B T 

T= 3hXW) ■ (8) 

This relation allows for an independent determination of the e-ph coupling strength 

A(c 2 ). 
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Figure 3. Evaluation of the temperature dependence of relaxation rates in EuFe2As2. 
(a) Time-dependent spectral weight of holes at fen > fcp for various temperatures. Data 
are vertically offset for clarity. Solid lines are exponential fits (see text), (b) Relaxation 
time constants Thoics determined from the data in (a) as function of temperature. The 
dashed line is a fit to equation (|8| in the temperature range T > 100K. 



Experimentally, the temperature-dependent relaxation of holes in the hole pocket 
near T in EuFe2As2 are used, which are found to be independent from the 
antiferromagnetic transition J26] . The respective hole relaxation time Thoies is determined 
from the temperature-dependent trARPES intensity by fitting exponential decay curves 
(fig. [3ta)) and is shown in fig. [3tb) as a function of temperature. The fit to equation (pi) 
for T > 100 K reveals a good agreement to the linear behavior and yields a value 

This 



of A (u 2 ) = 90(30) meV 2 , where the error represents a confidence interval of 

value, albeit considerably higher than the values found in the 3TM and in the evaluation 

of the excess energy, is in agreement with the other evaluations within error bars. 



3. Discussion 



The values of A (a; 2 ) determined by the three methods discussed above are compared in 
table [Tj Despite the limitations and drawbacks of the various methods, we find a perfect 
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Table 1. Values of A (w 2 ) determined by the three methods. 



compound T e (t) (A£(£)) r holes (t) 



EuFe 2 As 2 


56-65 


65(5) 


90(30) 


BaFei. 85 Coo.i5As 2 


46-55 


50(3) 


- 


BaFe 2 As 2 


30-46 


34(6) 


- 



agreement of all three methods within error bars. Especially, the trend for larger e-ph 
coupling in EuFe 2 As 2 than in BaFei. 8 5Co .i5As 2 and in BaFe 2 As 2 is nicely represented. 
Remarkably, all methods produce a small A {uj 2 ) < 100 meV 2 , which already indicates 
weak e-ph coupling. 

Our values of A {uj 2 ) can be compared to recently published values determined 
from optical pump-probe experiments. Mansart et al. [16] report a comparable value 
of A {uj 2 ) ~ 64meV 2 for Co-doped BaFe 2 As 2 , which was determined using the 3TM. 
Stojchevska et al. [25J derived a somewhat higher value of A {uj 2 ) = 110(10) meV 2 in 
SrFe 2 As 2 using temperature-dependent QP relaxation times, and for SmFeAsOi_. r F., ; 
an even larger value of A {uj 2 ) = 135(10) meV 2 is reported |27j . 

Based on our results of A {uj 2 ) we can estimate the value of the e-ph coupling 
constant A for a particular value of uj. Considering the Raman active Ai 9 mode at 
23 meV, which can be coherently excited [T7J |28j [29] and therefore shows enhanced 
e-ph coupling in the system, we find A < 0.2 for all compounds. This estimate is 
in agreement with calculations of various Fe-pnictide compounds, which report values 
of A < 0.35 [3U1 EI]. Taking the mean of the phonon spectrum as reference, A gets 
even smaller, in agreement with other publications [16, 25J. Even if we consider the 
lowest coupled modes around 12meV to be most important for e-ph coupling, A does 
not exceed a value of 0.5. Similarly small values for A have been found in the cuprate 
HTSCs [TU [15], which suggests limited importance of e-ph coupling for the pairing 
mechanism in both classes of materials. 

These boundaries for the e-ph coupling constant A allow to estimate the 
superconducting critical temperature T c , assuming a conventional BCS-type pairing 
based on e-ph coupling. Here, T c in isotropic systems and in a strong coupling regime 
is given by McMillan's formula [32J, modified by Allen and Dynes [33] . 



T-^expf 1-04(1 + A) ^ (9) 

Jc "l.20 eXP V A-/x*-0.62AW W 

where /i* is the effective Coulomb repulsion and uj\ og is the logarithmic average of the 
phonon spectrum. Taking /i* = and o;i g = 205 K [30], which provides an upper 
limit for T c , we find T c = 0.33K and T c = 7.5K for A = 0.2 and A = 0.5, respectively. 
Comparing these results with the experimentally found T c = 24 K in optimally doped 
BaFei.85Co .i5As 2 [31] and even up to T c = 38 K in K doped BaFe 2 As 2 [3S] demonstrates 
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that e-ph coupling in a conventional BCS pairing scenario cannot explain the high 
critical temperatures found in the pnictide HTSCs. In addition, albeit the stronger 
e-ph coupling, the critical temperatures found in EuFe2As2 upon Co [36J and K J37] 
substitution are considerably smaller compared to BaFe2As2. 

However, e-ph coupling might still play a significant role in the Cooper pair 
formation in the Fe pnictides. For instance, a very strong sensitivity of the Fe magnetic 
moment on the As height in the FeAs tetrahedra was found in DFT band structure 
calculations [3H], with a rate of 6.8 //b/A. Inelastic x-ray scattering experiments on 
CaFe 2 As2 concluded on strong coupling of phonons to magnetic excitations even in the 
high-temperature paramagnetic phase [2S]. Furthermore, a strong sensitivity of the 
maximum critical temperature on the pnictogen height was found pp. On the basis of 
these findings a strong magneto-phonon coupling in these compounds was proposed and, 
e.g. the strongly coupled K\ g mode perpendicular to the Fe layers, which modulates 
the pnictogen height, could mediate superconductivity in the Fe-pnictides [30]. Such a 
scenario is supported by the strong coupling of this particular mode to electronic states 
directly at the Fermi level [17] evidenced by the coherent excitation of this phonon 
mode [HI [283 EH]- 

4. Summary 

In summary, we demonstrated the quantitative analysis of the e-ph coupling strength 
from femtosecond time- and angle-resolved photoemission spectroscopy. In detail, we 
compared three different methods to determine the e-ph coupling strength in three 
122 FeAs compounds from time- and angle-resolved photoemission experiments. The 
transient temperature of the thermalized electronic distribution is analyzed in EuFe2As2, 
BaFei.85Coo.i5As2 and BaFe2As2 by a three-temperature model, while the rate of energy 
relaxation of the non-thermal electrons yields direct information on the second moment 
of the Eliashberg coupling function X(u 2 ). Finally, the temperature-dependence of 
hole relaxation rates in EuFe2As2 also allows to determine X(u 2 ). All three methods 
consistently yield a small A (u> 2 ) < 100 me V 2 , which results in A < 0.2 using a reasonable 
choice of the phonon spectrum. This value is discussed to be too small to explain 
superconductivity in the Fe pnictides by a conventional BCS-type pairing. 
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